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Abstract. We consider Lotka-Volterra systems in three dimensions depending on three real parameters. By 
using elementary algebraic methods we classify the Darboux polynomials (also known as second integrals) for 
such systems for various values of the parameters, and give the explicit form of the corresponding cofactors. 
More precisely, we show that a Darboux polynomial of degree greater than one is reducible. In fact, it is a 
product of linear Darboux polynomials and first integrals. 



1. Introduction 

The Lotka-Volterra model is a basic model of predator-prey interactions. The model was 
developed independently by Alfred Lotka (1925), and Vito Volterra (1926). It forms the basis 
for many models used today in the analysis of population dynamics. It has other applications 
in Physics, e.g. laser Physics, plasma Physics (as an approximation to the Vlasov-Poisson 
equation), and neural networks. In three dimensions it describes the dynamics of a biological 
system where three species interact. 

The most general form of Lotka-Volterra equations is 



oc £ £ i I ^^^^^ dir^jOCq^oCj j % 1^ 2 



We consider Lotka-Volterra equations without linear terms (e$ = 0), and where the matrix 
of interaction coefficients A = (a^) is skew-symmetric. This is a natural assumption related 
to the principle that crowding inhibits growth. 

The most famous special case of Lotka-Volterra system is the KM system (also known as 
the Volterra system) defined by 

(1) Xi = £j(xj+i — Xj_i) i = 1,2, ... ,n, 

where Xq = x n+ i = 0. It was first solved by Kac and van-Moerbeke in [13], using a discrete 
version of inverse scattering due to Flaschka [UJ. In [T7] Moser gave a solution of the 
system using the method of continued fractions, and in the process he constructed action- 
angle coordinates. Equations ([T]) can be considered as a finite-dimensional approximation 
of the Korteweg-de Vries (KdV) equation. The variables %i £11*6 cLll intermediate step in the 
construction of the action-angle variables for the Liouville model on the lattice. This system 
has a close connection with the Toda lattice, 



hi = ai(b i+1 - hi) i = 1, . . . , n - 1 
6j = 2(a--a-_ 1 ) i = l,...,n. 
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In fact, a transformation of Henon connects the two systems: 

1 , 

a% = -r\W2W i = l,...,n-l 

h = ~(x 2 i-i + x 2 i-2) i = l,...,n. 

The systems which we consider are all integrable in the sense of Liouville. In other words, 
there are enough integrals in involution to ensure the complete integrability of the system. 

Any constant value of a first integral defines a submanifold which is invariant under the 
flow of the Hamiltonian vector field. A second integral is a function which is constant on 
a specific level set. While a first integral satisfies / = 0, a second integral is characterized 
by the property / = Xf, for some function A which is called the cofactor of /. Second 
integrals are also called special functions, stationary solutions, and in the case of polynomials, 
eigenpolynomials, or, more frequently, Darboux polynomials. In systems which have a Lie 
theoretic origin (e.g. the full Kostant Toda lattice), they arise from semi-invariants of group 
actions. The importance of Darboux polynomials lies in the following simple fact. If / and 
g are relatively prime Darboux polynomials, with the same cofactor, then their quotient is 
a first integral. We propose to understand the behavior of a system based on the algebraic 
properties of its Darboux polynomials. 

As a starting point we consider the system 



ii = X\{rx 2 + sx 3 ) 

(2) x 2 = x 2 (-rx 1 +txs) 

£3 = %3(—sxi — tx 2 ) 

where r, s, t G R. 

Our main result is the following: 

Theorem 1. An arbitrary Darboux polynomial of the system is reducible. In fact, it is 
a product of linear Darboux polynomials. 

The method of proof that we use follows the approach of Labrunie in [T5] for the so called 
ABC system. 

The system fl2]) is Hamiltonian. We define the following quadratic Poisson bracket in M. 3 
by the formula 



. d d d d d d 

(3) 7T = rxxx 2 —— A h SX1X3— — A — \- tx 2 x 3 -— A — . 

OX\ OX 2 0X\ (JX3 0X 2 (JX3 

Generically, the rank of this Poisson bracket is 2 and it possesses a Casimir given by 
3. The function H = x\ + x 2 + X3 is always a constant of motion. In fact, taking 
H as the Hamiltonian and using the Poisson bracket ([3]) we obtain equations ([2]). 

Lotka-Volterra systems have been studied extensively, see e.g. [1], [13], [20]. The Darboux 
method of finding integrals of finite dimensional vector fields and especially for various types 
of Lotka-Volterra systems has been used by several authors, e.g. [6], [7], [15] [16], [18], [T9] . 

The paper is organized as follows. In Section 2, we recall a few basic facts about Darboux 
polynomials. In Section 3 we prove Theorem [1] under general conditions for r, s, t, and we 
also give the explicit form of the cofactors. Section 4 deals with the case s = t. We did 
not examine other such cases since the method of proof is identical with these two cases. 



Finally in Section 5 we present in detail three examples which include the open and periodic 
KM-system in three dimensions. 



2. Darboux polynomial preliminaries 
Consider a system of ordinary differential equations 

(4) ^ = Vi(xi(t), ...,x n (t)), i=l,...,n, 

where the functions Vi are smooth on a domain U C K n . Here K = M or IK = C, and we 
denote by K[x], x = (x\, ...,x n ), the ring of polynomials in n variables over K. Let <f> : I U 
be a solution of (jlj) defined on an open non-empty interval I of the real axis. A continuous 
function F : U — ► K. is called a /irsi integral of system (jlj) if it is constant along its solution, 
i.e. if the function F o is constant on its domain of definition for arbitrary solution <ft of 
(jlj). When F is differentiable, it is a first integral of system (jlj) if 

(5) L v (F) = 5>(x) — (x) = 0, 

i=l 4 

where L v is the Lie derivative along the vector field v = (vi, ...,v n ). If A is any function of 
x, then the Lie derivative of A is the time derivative of A, i.e. A = ^ = L V (A). The vector 
field generates a /?ou> 0t that maps a subset {7 of K n to K n in such a way that a point in U 
follows the solution of the differential equation. That is, 0(x)(t) = v(0(x)) V x G U. The 
time derivative is also called the derivative along the flow since it describes the variation of 
a function of x with respect to t as x evolves according to the differential system. 

Many first integral search techniques, such as the Prelle-Singer procedure, are based on 
the Darboux polynomials. A polynomial / G K[x] is called a Darboux polynomial of system 
© if 

(6) L v (f) = A/, 

for some polynomial A G K[x], which is called the cofactor of /. When A = 0, the Darboux 
polynomial is a first integral; / is said to be a proper Darboux polynomial if A 7^ 0. Let fx, 
/2 be Darboux polynomials with cofactors Ai, A2, respectively. It is easy to verify that: 

(i) The product /1/2 is also a Darboux polynomial, with cofactor Ai + A2, and 

(ii) If Ai = A2 = A then the sum fi + f 2 is also a Darboux polynomial, with cofactor A. 
The following propositions ([12]) give some more elementary but important properties of 
Darboux polynomials. 

Proposition 1. Let f,g& K[x] be non-zero and coprime (i.e. they do not have common 
divisors different from constants). Then, f\g is a rational first integral if and only if f and 
g are Darboux polynomials with the same cofactor A G K[x]. 

Proposition 2. (?) All irreducible factors of a Darboux polynomial are Darboux polynomials, 
(ii) Suppose that the system @ is homogeneous of degree m, i.e. all Vi are homogeneous of 
degree m, and let f be an arbitrary Darboux polynomial of @ with cofactor A. Then A is 
homogeneous of degree m — 1, and all homogeneous components of f are Darboux polynomials 
of (fjj) with cofactor A . 

Thus, the search for Darboux polynomials can be restricted to irreducible polynomials, 
and, if the system is homogeneous, to homogeneous polynomials. Since the dynamical system 



(j2J) is homogeneous of degree 2, the cof actor of any Darboux polynomial of the system will 
be a linear combination of the variables x±, x 2 , x 3 . It follows that any Darboux polynomial 
/ of system (j2J) will satisfy 

df df df 

(7) L(f) = x 1 (rx 2 + sx 3 )- ^x 2 (-rx 1 +tx 3 )- ^x 3 (-sx 1 -tx 2 )—- = (ax i +(3x 2 +'yx 3 )f, 

OX\ dx 2 dx 3 

where a, (3, 7 are constants. If it is not clear from the context, we shall denote these constants 
by a(f), (3(f), 7(/) respectively 

3. Darboux polynomials of the Lotka-Volterra system 

We carry out our analysis aiming at maximum generality, that is, imposing as few con- 
ditions on the parameters r, s, t as possible. In this section we make such assumptions in 
propositions [HI EH and Theorem [21 however, as we note in remark [31 one can obtain the 
results by making assumptions about the cofactor of the Darboux polynomial instead of the 
parameters. An important role in this work plays the homogeneity property, as can be seen 
in the following two propositions. 

Proposition 3. Let f be a homogeneous Darboux polynomial of system (Tj|) of degree m. If 
l{f) 0; then f has no x™ term so that f(x\,x 2 , x 3 ) = Xi(p(xi,x 2 , x 3 ) + x 2 ip(xi, x 2l x 3 ). 

Proof. Since the polynomial / is homogeneous, we use Euler's identity 

W df M df df 

(8) Xi- \-x 2 - \-x 3 -— = mf. 

OX\ dx 2 dx 3 

Using equation (jSJ) we substitute for in equation (CO) to obtain 

df df 

x 2 (-rxi+tx 3 -rx 2 -sx 3 )- \-x 3 (-sx 1 -tx 2 -rx 2 -sx 3 )- — = (axi + ((3-mr)x 2 +h-ms)x 3 )f. 

dx 2 dx 3 

Setting x\ = 0, x 2 = 0, and letting F(x 3 ) = /(0, 0, x 3 ) we have 

(9) - sx 2 3 F'(x 3 ) = (7 - ms)x 3 F(x 3 ). 

If s = 0, 7 7^ 0, equation (jHJ) implies that F — 0. Otherwise, if s 7^ we have F(x 3 ) = 
kx™ 7//s , for some constant k. Since / is homogeneous of degree m, the only term containing 
only x 3 is necessarily x™. Thus, if 7 7^ we must have F = also in this case, and the 
proposition is proved. □ 

We shall use the following notation: for a polynomial / = f(x±,x 2 ,x 3 ) we denote / = 
/L=o, /= f\x 2 =o, f = / |x3=o- We denote N m = {1,2, ...,m}, J\f m = N m U {0}, and for any 
number r, N m r = {nr : n G N m }. 

Proposition 4. Let f be a homogeneous Darboux polynomial of degree m. Ifj(f) 7^ 0, then 
s = =>- x 2 \f, t = 0=>x 1 \f, and 7 N m s x 2 \f, t 7^ and 7 N m t xi\f. 

We also have the following statements for a and (3: 

If (3(f) ^ 0, then r = x 3 \f, t = x^f, r ^ and (3 £ N m r x 3 \f, t ^ and 
(3<£-N m t^ Xl \f. 

If a (f) 7^ 0, then r = =>- x 3 \f, s = =>- x 2 \f, r ^ and a —N m r =>- x 3 \f, s 7^ and 
a ^ -N m s =>• x 2 \f. 



Proof. We prove the statements for 7. The proof of the statements for a and j3 is similar. 
Assume that 7 7^ 0. Then, it follows from Proposition [3] that / = x\4>i + X21P1, where 
0i = 4>i(%i,X2, x 3 ), tpi = ipi(xi,X2, x 3 ) are either homogeneous polynomials of degree m — 1, 
or zero (but they are not both zero). Setting this in equation (I7j) yields 

(10) xiL(0i)+x 2 L(^i) = (axi+f3x 2 + / yxs-rx2-sx 3 )x 1 (f) 1 + (ax 1 +f3x2+'yx3+rxi-txs)x2i>i. 
Setting x 2 = in equation (fit)]) we have 



x\L{4>\) = (axx + (7 - s)x 3 )x 1 (j)i. 

The operator 0i — > 0i commutes with the derivations with respect to X\ and x 3 , and therefore 
we obtain 

(11) «l^4) = («i+(7-.W*. 



sxix 3 \- — 

\OXi OX 3 

If s = then 0i = 0, which implies that 0i is divisible by X2 and that / = X\(j)\ + £2^1 is 
divisible by x 2 - Suppose now that s 7^ 0, deg0i = deg0i — m — 1, and that 7 7^ ns, n G iV m . 
The r.h.s. of (TTTj) is divisible by 27, and since 7 — s 7^ 0, it follows that Xi|0i. Let 0i = 2702, 
where 02 is a homogeneous polynomial of degree m — 2. Then, we have 

901 902 . 901 909 
h02, TT— =X r 



and from (TTTj) we obtain 

Since 7 — 2s 7^ 0, 02 is divisible by 27. Continuing in the same way we obtain 

/90 m _i 90 m _i\ 

\~9lc dx — / = ~ ~ l > s ) x VV m - x ' 

where deg0 m _i = 1, and Xi|0 m _i. Thus, m _i = const. Xi, and from the above equation we 
have sx 3 = ax\ + (7 — (m — 1)5)0:3. By equating coefficients we obtain 7 = ms, which is a 
contradiction. Therefore, we must have 0i = 0, which implies that / is divisible by x<i- 
Setting x\ = in (TTD1) and using (JTj) we obtain 

/ 9-01 dip 



tx 2 x 3 ^-^- - ^-i) = (13x2 + (7 - O^a)?! ■ 



If t = then ^1 = 0, hence -0i is divisible by 27 and so / is divisible by 27. Suppose that 
t 7^ 0, deg-01 = deg-01 = m — 1, and 7 7^ ni, n G ]V m . Then it can be shown in a similar 
way as above that -01 is divisible by 27, which implies that / is divisible by 27, and the 
proposition is proved. □ 

This leads to the characterization of the cofactors of Darboux polynomials of system (j2j), 
as follows. 

Proposition 5. Let f be a homogeneous Darboux polynomial of degree m. We have either 
lU) = 0, or 7(/) = 7is ; 71 G iV m; or 7(/) = 72*, 72 G iV TO; or 7(7) = 71s + 7^, 
72 G {1,2, ...,m - 1} ; 71 G iV m _ 72 . 



Proof. Since / is a Darboux polynomial it satisfies L(f) = (ax\ + (3x2 + 7^3)/- Suppose that 
7^0 and 7 7^ ns, n G iV m . Then by proposition @] / is divisible by X2, that is / = £2/1 for 
some homogeneous polynomial f\ of degree m—1 and we have 

L(h) = ((a + r)xi + (3x 2 + (7 - t)x 3 )fi. 

Suppose that 7(/i) 7^ 0, i.e. 7 7^ t, and that 7(/i) 7^ ns, n G iV m _i, that is 7 7^ ns + £, 
n G N m -\. Then, again by proposition 0] it follows that f\ is divisible by X2, and writing 
/1 = ^2/2 we obtain 

£(/ 2 ) = ((a + 2r)xi + /3a; 2 + (7 - 2t)x 3 )/ 2 . 

If 7 7^ 2£, and 7 7^ + 2t, n G iV m „ 2 , then / 2 is divisible by x 2 . Continuing in the same 
way, after m — 1 steps we obtain 

(12) £(/m-i) = ((" + (m - l)r)x! + #r 2 + (7 - (m - l)t)x 3 )/ TO -i, 

where deg/ m _i = 1. If 7 7^ (m — 1)£ and 7 7^ s + (m — l)t, then x 2 |/ m _i, and thus 
f m -i = const. x 2 . From equation (|T2*|) we then have — rxi + tx 3 = (a + (m — l)r)xi + /3x 2 + 
(7 — (m — l)t)x 3 , and by equating coefficients we obtain 7 = rat. We therefore conclude that 
we have either 7 = 0, or 7 = ns, or 7 = nt, n G iV m , or 7 = 71s + 7 2 t, 7 2 = 1, 2, m — 1, 
71 G N m - l2 , and the proposition is proved. □ 

Remark 1. We note that in the proof of proposition^ we can make the successive assump- 
tions 7(/) ^nt (ne iV„J, 7(/i) 7^ n£ ; (n G N m „i),..., 7(/ m -i) 7^ n*, (™ e fAicft imp/j/ 
i/iai £/ie respective functions are divisible by x\. We obtain the same result also in this case, 
in particular the relation 7is + 7 2 t with the conditions 71 = 1,2, ...,m — 1, 7 2 G A^ m _ 7l , which 
are the same with the conditions stated in the proposition. 

Proposition 6. Let f be a homogeneous Darboux polynomial of degree m. We have: 

(a) a(f) = 0, or a(f) = —air, cti G N m , or a(f) = —a 2 s, a 2 G N m , or a(f) = —a±r — a 2 s, 

a-2 = 1, 2, m - 1, a x G iV TO _a 2 . 

ft) /?(/) = 0, or /?(/) = ftr, A G iV W; or /?(/) = -fat, fa G jV W; or (3(f) = far - fat, 

fa = 1,2,..., 771-1, A G N m -p 2 . 

Proof. The proof is similar to the proof of proposition [5j 

The following propositions give further analysis on the cofactors, and their relation with 
the parameters and the form of the Darboux polynomials. 

Proposition 7. Let r, s, t be non-zero, r\s = q\, r\t = q 2 , and s\t = q 3 . Let f be a 

homogeneous Darboux polynomial of degree m, and a±, a 2; fa, fa, 71, 7 2 the integers which 
appear in propositions {5\ and 

(a) If ai + (a 2 -3)^ Af m -j and fa - (fa - j)± £ M m -j, for j = 0, 1, 2, m - 1, then 
ol 2 = fa. 

(b) If (a-i - j)q 1 + a 2 £ N m -j and 71 + (72 - j)± £ M m -j, for j = 0, 1, 2, m - 1, then 
"i = 72- 

(c) If -(fa - j)q 2 + fa i Afm-j and (71 - j)q 3 + 7 2 ^ N m -j, for j = 0, 1, 2, m - 1, then 
(3i = 7i- 

Proof. We prove statement (a) . The proof of statements (b) and (c) is similar. If a 2 or fa is 
non-zero, then by hypothesis we have 01(f) = — (a% + « 2 ^-) r 7^ and a(f) 7^ —nr, n G N m , 
or (3(f) = (fa — fa—)r 7^ and [3(f) 7^ nr, n G N m , respectively. In either case, it follows 



from proposition H] that / is divisible by X3. We can write / = £3/1, for some homogeneous 
polynomial f\ of degree m — 1, and we have 



L(/i) = ((a + s)xi + {(3 + t)x 2 + 7x3)/! 

= ((-air - (a 2 - l)s)xi + ((3 x r - (fo - l)t)x 2 + 7^3) /i- 

By the same argument as above, if we do not have 02(A) = /M/i) = 0> i-e. if we do not 
have «2 = /?2 = 1, then we have either a(fi) = —(at + (a 2 — l)^) r 7^ and a(fi) 7^ — nr, 
n G iV m _i, or = (/% - (/3 2 - l)^)r 7^ and /3(/ x ) ^ nr, n e iV m -i, and / x is divisible 

by X3. Continuing in the same way, after m — 1 steps we obtain 

(13) L(/ TO _i) = ((-air - (a 2 - (m - l))s)rr 1 + (ftr - (/3 2 - (m - l))t)ar 2 + 7x 3 )/m-i, 

where deg f m -i = 1- If we do not have a 2 = /3 2 = m — 1, then it follows by our assumptions 
that x 3 |/ m _i, which implies that / m _i = const. x 3 and / = const, x™. However, a 2 (x 3 ) = 
^2(^3) = 1, and by simple properties of Darboux polynomials it follows that a 2 (f) = P 2 (f) 
m. Therefore, we must have a 2 = (3 2 = n, for some integer n G {0, 1, 2, m}, and the 
proposition is proved. □ 

Proposition 8. Let r, s, t be non-zero, r\s = q±, and r\t = q 2 . Let f be a proper Darboux 
polynomial, homogeneous of degree m, with j(f) = 0, and let at, a 2 , f3i, j3 2 be the integers 
which appear in proposition® 

(a) If at\ 7^ and a.\q\ + a 2 ^ Af m , or (3\ 7^ and — (j3\q 2 — (3 2 ) ^ Af m , then s = —pt, for 
some positive rational number p. 

(b) If at = A = 0, (a 2 -j)± i iV m _. and -(/? 2 -j)i £ N m ^, for j = 0, 1, 2, m - 1, 
then we have f = x 3 2 I where I is a first integral. 

Proof, (a) Suppose a,\ 7^ and ctigi + a 2 ^ Af m . The other case is similar. Then we have 
a(f) = — (ai?i + a 2 )s 7^ and a(f) 7^ — ns, n G N m . From proposition H] it follows that / is 
divisible by x 2 , so that / = x 2 f\ for some homogeneous polynomial fi of degree m — 1, and 
we have 

(14) L(f 1 ) = ((a + r)x 1 +{3x 2 -tx 3 )f 1 . 

Equation ( IT^tl shows that f\ is a Darboux polynomial with 7(/i) = —t. However, from 
proposition [5] we have 7(/i) = 71 s + j 2 t for non-negative integers 71, 72 G {0, 1, 2, m — 1}. 
Therefore, 71S + 7 2 t = —t, which is possible only if ji 0, in which case s = _ an d 
the statement is proved with p = -±2£. 

(b) Suppose that a\ — fi\ — 0. Since / is a proper Darboux polynomial with "f(f) = we 
must have a 2 7^ or (5 2 7^ 0, and our assumptions imply that in fact a 2 = (3 2 (see proposition 
[7]). We have a(f) = —a 2 ^r 7^ and a(f) 7^ — nr, n G N m . It follows from proposition H] 
that / is divisible by £3. So / = x 3 f[ for some polynomial f[ of degree m — 1, and we have 

L{f[) = (-(« 2 - l)sx x - (ft - l)tx 2 )f[. 

By the same argument, if a 2 — 1 = j3 2 — 1 7^ 0, then is divisible by x 3 . Continuing in the 
same way, we find that / = x 3 2 I for some first integral I (I = 1 if a 2 = /3 2 = m), and the 
proposition is proved. □ 

These results allow us to characterize the Darboux polynomials of system (j2J). 



Proposition 9. Let f be a Darboux polynomial of system p| ; homogeneous of degree m. If 
s = then 

(15) f = xff 1} 

where fi is a Darboux polynomial with -y(fi) = 0. If s, t are non-zero and N m s n N m t = 0, 
then we have 

(16) / = xfxfh, 

where $2 is a Darboux polynomial with ^(f 2 ) = 0. Here, the non-negative integers 71, 72 are 
such that j(f) = 71s + 72^. 

Proof. If 7(7) = 0, then the result in each case follows by setting 71 = 72 = 0, fi — f 2 — f ■ 
Suppose that •y(f) 7^ and s = 0. Then, by proposition @] / is divisible by x 2 , and writing 
f = x 2 f[ we have 

L(f[) = ((a + r)x x + (3x 2 + (7 - t)x 3 )f[. 
Let this procedure be repeated as many times as it can, and let 72 be the number of times 
that it can. We have / = x 2 2 fi, where f\ is a Darboux polynomial with 7(/i) = 7 — 72^ = 
since we had to stop the division procedure by x 2 , and equation ( fl5|) is proved. Suppose 
now that ■y(f) 7^ 0, s, t are non-zero and N m s fl N m t = 0. Thus 7 ^ N m s or 7 ^ N m t. Let 
us consider the case 7 ^ N m s. The case 7 ^ N m t is similar. Then / is divisible by x 2 and 
as before we have / = x 2 2 f 2 , where f 2 is a Darboux polynomial with , ~f(f 2 ) = 7 — 72^- Since 
we had to stop the division procedure by x 2 , we must have either i(f 2 ) = 0, in which case 
equation (jToT) is satisfied with 71 = and f 2 = f 2 , or j(f 2 ) — 7i s > f° r some 71 G N m . In the 
latter case j(f 2 ) ^ N m t and f 2 is divisible 71 times by xi, that is, f 2 = xj 1 f 2 , i(f 2 ) = 0, and 
equation ({ToT) follows. □ 

Remark 2. T/ie condition N m s fl iV m £ = implies that there do not exist integers tii, 
n 2 G N m such that s = —t. This condition is satisfied in each of the following cases: 

(a) one of s, t is positive and the other is negative, 

(b) s, t have the same sign but one is rational and the other irrational, 

(c) s, t have the same sign, they are both irrational, and their ratio is irrational, 

(d) s, t have the same sign, they are both rational, and s/t < 1/m or s/t > m, 

(e) s, t have the same sign, they are both irrational, their ratio is rational, and s/t < 1/m 
or s/t > m. 

Remark 3. In proposition instead of the condition N m s fl N m t = 0, we can make an 
alternative assumption as follows. First let s\t = q 3 , and let f k , k = 0, 1, 2, / = f , be a 
sequence of Darboux polynomials as we describe below. We denote i(f k ) = li(fk)s + i 2 (fk)t, 
7i = 7i(/)> 72 = 72(/)- For k = 0, 1,2, ...,71 + 72 - 1, we suppose that 

(17) (i) 7i(/ fc ) + j 2 (f k ) — £ Nm-k or (ii) ji(f k )q 3 + 72(A) £ Kn-k- 

In particular, if^i(f k ) = then we require condition (i) to hold, whereas if^ 2 (f k ) — then we 
require condition (ii) to hold [if'ji(fk) 7^ and^ 2 (f k ) 7^ then we can have either condition 
(i) or (ii)). If condition (i) holds, then j(f k ) 7^ and j(fk) 7^ ns, n G N m - k , which implies 
that f k is divisible by x 2 . Thus f k = x 2 f k+1 , and ji(f k +i) = li(fk), J2(f k +l) = 12(h) - 1- 
// condition (ii) holds, then j(f k ) 7^ and j(f k ) 7^ nt, n G N m ^ k , which implies that X\\f k . 
In this case we have f k = Xxf k+1 , ji(f k +i) = Ji(f k ) ~ 1, ^(fk+i) = 72 (fk)- Following this 
procedure, after 71 + 72 steps we obtain f = x^x], 2 f , where i(f') = 0. 



The following proposition states similar results in terms of the constants a and fi. The 
proof is similar to the proof of proposition [9j 

Proposition 10. Let f be a homogeneous Darboux polynomial of degree m, and let ax, 

fix, @2 be the integers which appear in proposition® 

(i) If s = then f = x^fx, where fx is a Darboux polynomial with a(fx) = 0. 

(ii) Ifr, s are non-zero and —N m rn (—N m s) = then f = x^x^ 2 f 2 , where / 2 is a Darboux 
polynomial with a(f 2 ) = 0. 

(Hi) If r, t are non-zero and N m r R (—N m t) = then f = x^x^ 2 fs, where / 3 is a Darboux 
polynomial with fi{fz) — 0. 

We are now ready to prove the main result of this section. 

Theorem 2. Let f be a Darboux polynomial of system (01) ; homogeneous of degree m. 
Suppose that either: (i) s = and N m r n (—N m t) = ; or (ii) r, s, t are non-zero, 
N m r n (-N m t) = 0, N m s n N m t = 0, and (-N m r) H (-N m s) = 0. (In particular, con- 
dition (ii) is satisfied, for example, when r > 0, t > and s < 0, orr<0,t<0 and s > 0). 
Then, there exist three non-negative integers i, j, k and a polynomial first integral I -which 
may be trivial- such that 

(18) / = x[xix k 3 I 
and 

(19) a(f) = -rj-sh, fi(f) = ri-tk, -y(f) = 8i + tj. 

Proof. Consider the case s = and N m r D (—N m t) = 0. The other case is similar. We use 
equation (fTS]) of proposition and the equations in statements (i) and (iii) of proposition [TU] 
in the following algorithm. 

(1) Set n = and /„ = /. 

(2) Applying proposition [TUl for a (statement (i)) yields 

fn = xTfn+l, a(f n+1 )=0. 

If f n+ x is a first integral, go to the final step, else increment n by one. 

(3) Applying proposition [TUl for fi (statement (iii)) yields 

fn = X^X^Ux, fi(f n+ x) = 0. 

If f n+ \ is a first integral, go to the final step, else increment n by one. 

(4) Applying proposition [9] for 7 (eq. (TT5l) ) yields 

fn = Xffn+x, 7(/rH-l) = °- 

If f n+ x is a first integral, go to the final step, else increment n by one and return to step 2. 

(5) (Final step) Set I = f n+1 and using the sequence of equations linking fi to 1 = 1, n 
given by the algorithm determine the exponents i, j, k in eq. (|18|) . 

At every step one has deg/z + i < deg/;; when three consecutive terms of the sequence are 
of the same degree, they are equal and a(fi) = fi(fi) = j{fi) = 0, so fi is a first integral. 
Thus the algorithm converges in a finite number of steps. Equation ( jT9l) follows from simple 
properties of Darboux polynomials. 

If condition (ii) holds, then the proof is the same but now in steps 2 and 4 of the algorithm 
we use the equation in statement (ii) of proposition [10], and equation ( |T6l) of proposition [9J 
respectively. □ 



4. The case s = t 



In this section we study the case s = t, which is not covered by Theorem [2] in the previous 
section. It can be seen that in this case x\ + is an additional linear Darboux polynomial 
of system ([2]), with cofactor sx^. Therefore, polynomials of the form / = x\x\x\[x,\ + X2) 1 , 
where i, j, k, I are non- negative integers, are Darboux polynomials. We show that a Darboux 
polynomial will have this form with I > 0, provided its cofactor satisfies some conditions 
which depend on the ratio r\s. 

Proposition 11. Suppose that r, s, t are non-zero, s = t, and let r\s = q\. Let f be a 

homogeneous Darboux polynomial of degree m which does not have the form (1 1 S\) . and let a±, 
a 2, Pi, Pi, 7i; 72 be the integers which appear in propositions^ and\^ For j = 0, 1, 2, m — 
1, suppose that a x + (a 2 - j)± £ N m -j, Pi - (P2 ~ j)j^ £ N m -j, ("i - j)Qi + "2 i Nm-j, 
and —((Pi — j)qi — P2) ^ Nm-j- Then, we have (i) a 2 = P2 and (ii) ai + Pi < 71 + 72- 

Proof. Relation (i) is statement (a) of proposition [71 We prove the inequality (ii). Suppose 
on the contrary that a± + Pi > 71 + 7 2 . By arguments that we have used repeatedly in this 
paper (for example see proposition [H]), / is divisible ai times by x 2 and Pi times by X\. Thus 
we have / = x^x^ 1 f, where /' is a Darboux polynomial of degree m — (ati + Pi) such that 

L (f) = {-®2SXi - p 2 sx 2 + (7 - ("i + Pi)s)x 3 )f. 

By proposition [5] there exist non-negative integers j[, 72 G {0,1,. ..,m — (ai + Pi)} such 
that j(f') = JiS + 7 2 t = (j'i + 7 2 )s. This implies that j[ + 7 2 = 71 + 72 — a.\ — Pi < 0, a 
contradiction. 

If «i + Pi = 71 + 72, then from the equation above we have L(f') = (— a 2 sxi — P2sx 2 )f, 
and our assumptions imply that /' is divisible a 2 times by x 3 (proposition [H]) . So we have 
/' = X3 2 /, and therefore / = J, where / is a first integral. Since we assume that / 

does not have the form ( 1181) we may exclude this possibility, and the proof is completed. □ 

Proposition 12. Suppose that r, s, t are non-zero, s = t, and let r\s = q\. Let f be a 

homogeneous Darboux polynomial of degree m which does not have the form (\18§ . With the 
same assumptions as in proposition 1771 we have f = (xi + a^/i, for some polynomial fi. 

Proof. From proposition [8] it follows that we may assume 7(/) 7^ 0. By proposition [3] / does 
not have an x™ term and we can write / = x\<§>\ + x^i, for some polynomials <pi, ipi. For 
a polynomial / = f(xi,X2,x%) we denote by / the polynomial obtained from / by setting 
X2 = ~xi z thatjs / = f(xi,x 3 ) = f\ X2 =- xl . So, / = x t (^i - ipi), and letting h x = <fa.- V»i 
we have / = xihi. Setting s = t and X2 = —x\ in equation ( TTDT) we obtain 

x x L((j) X ) - xiL(ipi) = ((a- p + r)xi + (7 - s)x 3 )xi($ 1 - $1) 

or 

(20) L^) = ((a-p + r)xi + ( 1 -s)x z )hi. 

Setting s = t and X2 = —x\ in equation (JTj) we obtain 



(21) 



Combining equations ( !20l) and ( I2T1) . and noting that a 2 = ft (proposition [TTJ), we obtain 

,/dhi dhi\ 

(22) - xi(rxi - sx 3 ) - —J = (-(oil + ft - l)rx x + (71 + 72 - l)sx 3 )%- 

From proposition [TT] we also have a>i + ft < 71 + 72, which implies that the term — (cti + 
ft — l)rxi + (71 + 72 — l)sx 3 is not a constant multiple of (rxi — sx 3 ). Since (rxi — sx 3 ) 
divides the right-hand side of equation (1221) . it divides h\. Therefore we have 

hi = (rxt - sx 3 )p 1 + (-rx 2 - sx 3 )xi, 

for some polynomials p%, Xi- Let h 2 = Pi + Xi- Then, hi = (rxi — sx 3 )h 2 and / = 
Xi(rxi — sx 3 )h 2 . We have 

(23) 7— = (rxi - sx 3 )- h rpi + (rxi - 5x3)7: — , 

axi axi axi 

/ 0/1 x d^i / .dpi . dxi 

(24) —— = [rxi - sx 3 )- V{rxi-sx 3 )- r\i- 

ox 2 dx 2 dx 2 

Substituting for from equations ( 1231) . ( 1241) respectively in equation ( 1221 we obtain 

( \(t \( d P l , d * l \ t \( d P 1 , d * l \ , 

-X!(rx! - ,x 3 )((rx 1 - sx 3 ) + ^-J - (ra* - sx 3 )[— + — ) + 

+r(pi + xi)) = (-(«i + fh - l)rxi + (71 + 72 - 1)^3)^1 
and simplifying further we have 

(25) - xi(rxi - sx 3 ) - -^-pj = (-(ai + ft - 2)rxi + (71 + 72 - l)sx 3 )A 2 . 

The term — («i + ft — 2)rxi + (71 + 72 — l)sx 3 is not a constant multiple of (rxi — sx 3 ), and 
so (rxi — 5x3)1/12- Continuing in the same way we find that / is divisible by an infinity of 
powers of (rxi — SX3) , which is a contradiction. Therefore we must have / = 0. This implies 
that / = (xi +x 2 )fi, for some polynomial fi, and the proof of the proposition is completed. 
□ 

Corollary 1. Suppose that r, s, t are non-zero, s = t, and let r\s = q x . Let f be a 
homogeneous Darboux polynomial of degree m. With the same assumptions as in proposition 
[771 we have 

(26) f = x\xlx k 3 (xi + x 2 ) l I, 

where I is a first integral and i, j , k, I are non-negative integers. 

Proof. Note first that if 7(/) = then by proposition [8] it follows that we must have a,\ = 
ft = and / = X3 2 /. If 7(/) 7^ and / does not have the form (JTBJ) (which is (I2"B"]) with 
/ = 0), then by proposition [T2] we have / = (xi + x 2 )/i for some polynomial fi, and 

L(/i) = (axi + ftr 2 + (7 - s)x 3 )fi. 

Repeating this procedure a finite number of steps, we find that / has the form f[26|) . □ 



Remark 4. Similar results hold when r = s and r = —t. It can be seen that if r = s then 
£2 + 2:3 is a linear Darboux polynomial with coj "actor —rx\. Under conditions analogous to the 



ones we have used in this section, we have f 



x\x 2 X 3 { 



X2+X3) I. Similarly, if r 



-t then 



X1+X3 is a linear Darboux polynomial with cof actor rx 2 , and we have f = x\x 2 x^(xi + x 3 ) k I. 



5. Examples 

5.1. Example: The KM-system. We give a complete description of Darboux polynomials 
for the case of the KM system (r = 1, s = 0, t = 1). 



(27) 



X\ = Xix 2 

X 2 = —X\Xi + x 2 x 3 

x 3 = -x 2 x 3 



The Hamiltonian description of system (QQ) can be found in [10] and [8] . We will follow 
and use the Lax pair of that reference. The Lax pair in the case n = 3 is given by 



L=[B,L), 



where 





Xi + x 2 
■sJX\X 2 







V 







y/X 2 X 3 



yJX Y X 2 \ 

y/X 2 X 3 
X 2 + X 3 

x 3 J 



and 



B 



\ 






\^x x x 2 

~\y/x 2 x 3 



\y/X X X 2 \ 

\^x 2 x 3 

" 

/ 



This is an example of an isospectral deformation; the entries of L vary over time but the 
eigenvalues remain constant. It follows that the functions Hi = }TrL* are constants of 
motion. We note that 

Hi = 2(xi + x 2 + x 3 ) 
corresponds to the total momentum and 

3 2 

h 2 = ^2 x i + 2 Xi%i + i ■ 



Using (jSJ) we define the following quadratic Poisson bracket, {xi,x i+ i} = XiX i+ i,i = 1,2, 
and {xi, x 3 } = 0. For this bracket detL = Casimir and the eigenvalues of L are in 

involution. Taking the function Hi = X\ + x 2 + £3 as the Hamiltonian we obtain equations 
( 1271) . Therefore the system has a Casimir given by F = xix 3 and a constant of motion 
a^i + £2 + x 3 corresponding to the Hamiltonian. Note that H 2 = Hf — 2F. 

In the following table we present all Darboux polynomials of degree < 3 and the corre- 
sponding cofactors. 





Darboux polynomial 


cofactor 


-1 

1 


X\ 


X2 


2 


X 2 


-Xl + x 3 


o 




X 3 


-x 2 


A 

4 


Xi+ x 2 + x 3 


U 



Table 1. Linear Darboux polynomials and corresponding cofactors 





Darboux polynomial 


cofactor 




Darboux polynomial 


cofactor 


1 


rp ^ 


2x 2 


6 


x 2 x 3 


-Xi - x 2 + x 3 


2 


rp ^ 


— 2xi + 2x 3 


7 


xi(xi + x 2 + x 3 ) 


x 2 


3 


2 
x 3 


-x 2 


8 


x 2 (xi + x 2 + x 3 ) 


-Xi + x 3 


4 


X±X2 


-Xi + x 2 + x 3 


9 


x 3 (xi +x 2 + x 3 ) 


-x 2 


5 


XiX 3 





10 


Ci(x{ + x 2 + x 3 + 
+2xix 2 + 2x 2 x 3 ) + c 2 x\x 3 






Table 2. Quadratic Darboux polynomials and corresponding cofactors. Note that (10) is 
a sum of two first integrals, and thus a first integral; c\, ci are constants 





Darboux polynomial 


cofactor 




Darboux polynomial 


cofactor 


1 




3x 2 


9 


X^Xi + X 2 + X 3 ) 


— 2xi + 2x 3 


2 


x 2 


— 3xi + 3x3 


10 


X^(xi + X 2 + X 3 ) 


-2x 2 


3 


3 
X 3 


-3x 2 


11 


XiX 2 (xi + X 2 + X 3 ) 


-Xi + x 2 + x 3 


4 


x\x 2 


-Xi + 2x 2 + x 3 


12 


^2^3(^1 + x 2 + X 3 ) 


-Xi - x 2 + x 3 


5 


x\xi 


-2xi + x 2 + 2x 3 


13 


CiXi(x^ + x\ + X3 + 

+2xix 2 ) + C2X^x 3 


x 2 


6 


x l 2 x 3 


-2xi - x 2 + 2x 3 


14 


C 3 X 2 (x'f + x\ + X3 + 

+2x x x 2 + 2x 2 x 3 ) + C4X1X2X3 


-Xi + x 3 












7 


x{x 2 


-xi - 2x 2 + x 3 


15 


£5^3(^1 + ^2 + ^3 

+2xix 2 + 2x 2 x 3 ) + 06X3X1 


-x 2 


8 


x\{x\ + x 2 + x 3 ) 


2x 2 









Table 3. Cubic Darboux polynomials and corresponding cofactors; c±, ...,Cq are constants 



5.2. Periodic KM-system. The periodic system (r = 1, s — —1, t — 1) 

The periodic KM-system is given with the same equations ([T|) plus a periodicity condition 
Xi = Xi +n . In the case n = 3 we obtain: 

Xi = X1X2 — X1X3 

x 2 = -XiX 2 + x 2 x 3 
X 3 = X1X3 - X 2 X 3 

We give a different type of Lax pair for this system from pQ. 

/0 xi 1\ 
L = 1 x 2 
1 / 



(0 £l£2\ 

x 2 x 3 . 
X3X1 / 

It follows that the functions Hi = }TrL* are constants of motion. We note that H\ = , 
H 2 = £1 + £2 + X3 and i?3 = 1 + £1X22:3. As expected the function H 2 = x\ + £2 + £3 plays 
the role of the Hamiltonian with respect to the poisson bracket ([3]) while F = X1X2X3 is a 
Casimir. 

In the following table we present all Darboux polynomials of degree < 3 and the corre- 
sponding cofactors. 





Darboux polynomial 


cofactor 


1 


£1 


x 2 - x 3 


2 


x 2 


-X X + X 3 


3 


x 3 


Xi - x 2 


4 


Xl + £2 + x 3 






Table 4. Linear Darboux polynomials and corresponding cofactors 





Darboux polynomial 


cofactor 




Darboux polynomial 


cofactor 


1 


,->-» ^ 


2x 2 — 2x 3 


6 


x 2 x 3 


-x 2 + x 3 


2 


,->-» ^- 

x 2 


-2xi + 2x 3 


7 


Xi(Xi + x 2 + £3) 


x 2 - x 3 


3 


•> 

™^ 

x 3 


2xi — 2x2 


8 


x 2 (xi + x 2 + £3) 


-Xi + x 3 


4 


X1X2 


—X\ + x 2 


9 


x 3 (xi + x 2 + £3) 


Xi — X2 


5 


X1X3 


Xl - x 3 


10 


(xi + x 2 + x 3 ) 2 






Table 5. Quadratic Darboux polynomials and corresponding cofactors 





Darboux polynomial 


cofactor 




Darboux polynomial 


cofactor 


1 


x\ 


3x 2 — 3x 3 


11 


x 2 (xi + x 2 + x 3 ) 2 


-Xi + x 3 


2 


A 


— 3xi + 3x 3 


12 


x 3 (xi +x 2 + x 3 ) 2 


Xi - x 2 


3 


x\ 


3xi — 3x 2 


13 


x'{(xi + X 2 + £3) 


2x 2 - 2x 3 


4 


x\x 2 


-xi + 2x 2 - x 3 


14 


x\{xi + X 2 + £3) 


-2xi + 2£ 3 


5 


X1X3 


xi + x 2 - 2x 3 


15 


x'l(xi + X 2 + £3) 


2xi — 2x2 


6 




— 2xi + x 2 + x 3 


16 


XiX 2 (Xi + £ 2 + x 3 ) 


-Xi + x 2 


7 


2 

X 2 X3 


-xi - x 2 + 2x 3 


17 


x 1 x 3 (x 1 + £ 2 + x 3 ) 


Xi — x 3 


8 


2 

X3X1 


2xi — X2 — X3 


18 


x 2 x 3 (xi + x 2 + £3) 


-£ 2 + x 3 


9 


X3X2 


xi - 2x 2 + x 3 


19 


£l£ 2 £ 3 





10 


xi(xi + x 2 + x 3 y 2 


x 2 - x 3 


20 


c\{x\ + £2 + X3 + 3x 2 x 2 + 

+3x^X3 + 3X2Xi + 3X2X3 + 

+3x|xi + 3x|x2) + C2X1X2X3 






Table 6. Cubic Darboux polynomials and corresponding cofactors; c\, C2 are constants 



5.3. The case s = t (r = 5, s = t = 1). 

i\ = 5xiX 2 + Xix 3 

x 2 = -5xix 2 + x 2 x 3 
X 3 = -X1X3 - x 2 x 3 

We list below all linear, quadratic, and cubic Darboux polynomials of the above system 
which do not have the form (ITS]) , and their corresponding cofactors. 





Darboux polynomial 


cofactor 


1 


Xi + x 2 


£3 



Table 7. Linear Darboux polynomials and corresponding cofactors 





Darboux polynomial 


cofactor 




Darboux polynomial 


cofactor 


1 


Xi(xi + x 2 ) 


52 2 + 22 3 


4 


(21 + 2 2 ) 2 


22 3 


2 


x 2 (x 1 + x 2 ) 


-52i + 22 3 


5 


(2i + 2 2 )(2i + 2 2 + 2 3 ) 


23 


3 


£3(21 + 2 2 ) 


-2i - X 2 + 2 3 









Table 8. Quadratic Darboux polynomials and corresponding cofactors 





Darboux pol. 


cofactor 




Darboux polynomial 


cofactor 


1 


(21 + 2 2 ) S 


32 3 


9 


2i2 3 (2i + 2 2 ) 


-2i + 42 2 + 22 3 


2 


2i(2i + 2 2 ) 2 


52 2 + 32 3 


10 


2 2 2 3 (2i + 2 2 ) 


-621 - 2 2 + 22 3 


3 


22(21 + 2 2 ) 2 


— 52 1 + 323 


11 


2i(2i + 2 2 )(2i + 2 2 + 2 3 ) 


52 2 + 22 3 


4 


23(21 + 2 2 ) 2 


—X\ - X 2 + 22 3 


12 


2 2 (2i + 2 2 )(2i + 2 2 + 2 3 ) 


— 52i + 223 


5 


2 2 (2i + 2 2 ) 


102 2 + 323 


13 


23(21 + 2 2 )(2i + 2 2 + 2 3 ) 


— 2i - 2 2 + 2 3 


6 


23(21 + 2 2 ) 


— IO21 + 32 3 


14 


(2i + 2 2 )(2i + 2 2 + X 3 ) 2 


x 3 


7 


X 2 3 {Xi +2 2 ) 


-22i - 22 2 + 2 3 


15 


(2i + 2 2 ) 2 (2i + 2 2 + 2 3 ) 


22 3 


8 


2i2 2 (2i + 2 2 ) 


— 52 1 + 52 2 + 323 









Table 9. Cubic Darboux polynomials and corresponding cofactors 
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